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Abstract. We investigate cosmologies where the accelerated expansion of the Universe 
is driven by a field with an anisotropic equation of state. We model such scenarios within 
the Bianchi I framework, introducing two skewness parameters to quantify the deviation 
of pressure from isotropy. We study the dynamics of the background expansion in these 
models. A special case of anisotropic cosmological constant is analyzed in detail. The 
anisotropic expansion is then confronted with the redshift and angular distribution of 
the supernovae type la. In addition, we investigate the effects on the cosmic microwave 
background (CMB) anisotropics for which the main signature appears to be a quadrupole 
contribution. We find that the two skewness parameters can be very well constrained. 
Tightest bounds are imposed by the CMB quadrupole, but there are anisotropic models 
which avoid this bound completely. Within these bounds, the anisotropy can be beneficial 
as a potential explanation of various anomalous cosmological observations, especially in 
the CMB at the largest angles. We also consider the dynamics of linear perturbations in 
these models. The covariant approach is used to derive the general evolution equations 
for cosmological perturbations taking into account imperfect sources in an anisotropic 
background. The implications for the galaxy formation are then studied. These results 
might help to make contact between the observed anomalies in CMB and large scale 
structure and fundamental theories exhibiting Lorentz violation. 
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1. Introduction 

There is a remarkable amount of observational evidence that the large-scale structure of the 
Universe resembles nearly the simplest and most symmetric imaginable system. Hence, it 
is well approximated by the flat geometry version of the Friedmann-Lemaitre models which 
are both homogeneous and isotropic. The fluctuations in the cosmic microwave background 
(CMB) are a crucial probe of these properties, and therefore the recent detections of some 
unexpected features in the CMB temperature anisotropics have raised a lot of speculations 
about the need to reconsider some of the basic cosmological assumptions. A hemispherical 
asymmetry has been reported [1]. The angular correlation spectrum seems to be lacking 
power at the largest scales [2] . The ahgnment of the quadrupole and octupole (the so called 
Axis of Evil [3]) could also be seem as an extra-ordinary and unlikely result of statistically 
isotropic perturbations, even without taking into account that these multipoles happen 
also to be aligned to some extent with the dipole and with the equinox. The Axis of Evil 
does not show any correlation with the lack of angular power [4]. Given the a posteriori 
nature of these considerations, the cosmic variance (and the inevitable arbitrariness in any 
statistics), the statistical significance of all these anomalies is indeed a debatable issue 
[5-8]. 

Various possible cosmological effects have been proposed as possible explanations for 
those anomalous features. For instance, to introduce a preferred axis, one has to generate 
some isotropy breaking. According to whether this occurs at an early time or at late times, 
one may classify the models into those in which an originally isotropic CMB fluctuation 
is distorted on its way to us [9] and into those whose statistical anisotropy is imprinted 
in the primordial fluctuations [10]. A primordial imprint has been suggested to be left by 
isotropization happening during inflation [11, 12], parity violating couplings [13] vector field 
gradients [14] or non-standard spinors [15]. On the other hand, the apparent alignments 
could also be caused by local effects in our neighborhood. Notice, however, that the Axis 
of Evil has been argued to probably persist after foreground removal [16]. Local voids have 
been though suggested as its possible origin [17] (and quite interestingly such voids could 
explain the present days acceleration too [18]). 

In the case CMB is distorted during the late acceleration, the signatures of anisotropy 
would automatically be seen at the smallest multipoles of the CMB, since the perturbation 
wavelengths corresponding to these angles enter inside the horizon at the same epoch that 
the dark energy dominance begins. Therefore one would have one coincidence problem 
less. The apparent statistical anisotropy has been associated with dark energy in the form 
of vector fields [19-22], shear viscous fiuids [23] and elastic sohds [24]. See also [25-28] 
for other investigations. In general, if dark energy is something else than the cosmological 
constant A-term, one expects it to have anisotropic stresses at least at the perturbative 
level. This is also a generic prediction of modified gravity theories [29-35]. Moreover, it 
could be used to distinguish among scalar field theories in which the possible nonminimal 
coupling occurs only at the perfect fluid matter sector [36-39]. Such perturbative shear 
stresses have in general been attempted to be constrained through parameters describing 
directly the difference of the metric variables [40] and the viscous properties of generalized 
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fluids [41-43]. Unfortunately, the perturbative anisotropic stress could escape detection, 
when its influences are restricted to so large scales that they are only seen through the 
amplitude of the first few multipoles of the CMB. In spite of that, in many specific models 
the Jeans scale of the (effective) dark energy is small enough for dark energy to form small- 
wavelength perturbations which might be probed by e.g. weak lensing experiments [44- 
47]. Moreover, a statistical anisotropy due to dark energy would leave specific signatures 
which are not described alone by the amplitude of the CMB angular power spectrum 
~ {aim, ci'^ra) ^ whcrc aim are the coefficients of the spherical expansion of the anisotropy. 
Instead, correlations between different multipoles / are predicted, and also m-dependent 
patterns (whereas statistical isotropy would predict, on average, the same value regardless 
of m). The theoretical prediction from statistical anisotropy to the CMB is thus the 
nondiagonality of the matrix {aim, "^Li)- 

Here we continue a previous study which investigated the implications and origins 
of an anisotropic acceleration [48]. There we studied exact anisotropic but homogeneous 
solutions of the cosmological equations. We now extend such investigation and generalise 
to the case where inhomogeneities are present and introduce a covariant framework for 
the study. The Bianchi I model is sufficiently simple to allow semi-analytical calculations, 
and it captures the basic effect present also in more complicated anisotropic models by 
featuring their important common property, direction-dependent expansion rates. Such 
metric is useful to describe magnetic fields [49] or their Yang- Mills generalizations [50] . As 
is well known, the anisotropy due to such fields tends to decay. Usually the constraints 
inferred on the Bianchi type I models have indeed been for cases in which the universe may 
initially be anisotropic but will then isotropize. Solutions which do not isotropize would 
require anisotropic matter sources, which in addition should dilute slower than radiation 
or dust. There has been no evidence that such cosmological sources would exist. However, 
as the universe is presently believed to be dominated by some source with unexplained 
negative pressure, a question we would like to ask is whether this pressure could also be 
anisotropic. This is of interest even without the apparent hints of statistical isotropy at 
the CMB. 

After the study of Barrow on constraining anisotropic stresses in the late universe [49] , 
we note that very recently there has been also interest in employing Bianchi I to study the 
CMB. CampaneUi et al [51] describe magnetic fields present at the last scattering epoch 
resulting in an ellipsoidal universe. Similar effects have been proposed due to moving 
dark energy [52]. Rodrigues has suggested that an anisotropic cosmological constant 
[53] could result from an infrared noncommutative property of the spacetime. Longo 
has claimed that handedness of spiral galaxies is anisotropic [54] and that this, (together 
with the Axis of Evil) could be a result of large-scale magnetic fields [55]. If true, this 
would be a strong further motivation for our study since our model could unify magnetic 
and dark energy fields by describing Yang-Mills field or the mutually coupled system. 
In addition to these approaches, the anisotropic expansion rates have been exploited in 
infiationary considerations [11, 12, 15, 56, 57]. Then the appearance of anomalous features 
can be consistent, but does not necessitate, the existence of anisotropic sources. Infiation 
can smooth out initial anisotropies and inhomogeneities, but if the universe did not undergo 
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too many e- folds of inflation [58,59], some traces of such possible anisotropics could now 
be detectable at the largest scales. On the other hand, if one assumes isotropic initial 
conditions for the inflation, anisotropies could originate during inflation from a nonstandard 
energy source. This could be the same as dark energy. 

We attack the problem at three fronts: phenomenologically, theoretically and 
observationally. In section 2 we present our phenomenological description of anisotropic 
sources in terms of the skewness parameters 5 and 7, and derive the basic equations 
describing the model as a dynamical system. We then check the generic asymptotic 
evolution of the universe in some simple cases. This is a unifying description which can 
be used for any particular model of inflation or dark energy. Especially we contemplate 
the possibility of generalizing the usual A term in such a way that it would correspond to 
a constant vacuum energy that may exert anisotropic pressure. Possibilities of unifying 
inflation and dark energy arise therefore naturally. Section 4 is devoted to the study of 
the implications of these models to cosmological observations. Explicit constraints are 
derived from the luminosity distance-redshift relationship of the supernovae of type la 
(SNIa) and the amplitude of the quadrupole anisotropy in the CMB. We find that the 
additional parameters we introduced, 5 and 7, are well constrained. Finally, the last part 
of this article considers the inhomogencitics in the models considered here. In Section 5.1 
wc briefly introduce our notation for the covariant formalism and then present the general 
equations describing anisotropic cosmologies. Then we apply these general results into 
some interesting cases. In section 5.2 we study then large scale structure formation under 
various specific assumptions about the model. We conclude in section 6. 

2. Anisotropic Equation of State 

We set up our framework within the Bianchi type I case. Usually the Bianchi models 
isotropizc [60]. This is, in particular, true of the Bianchi I model with a perfect fluid 
content. Models that do not isotropize, have not been studied to such extent (though, 
for instance Ref. [61] presents exact solutions of inflationary Bianchi I cosmologies in the 
presence of various anisotropic sources). There has been good reasons for that: firstly, 
the modern cosmological observations, most importantly the satellite measurements of the 
CMB have established that the universe at large scales is isotropic to about one part in 
10^. Secondly, one does not usually expect growing anisotropies, on the contrary to have 
them usually requires some imperfect and thus exotic matter sources. Thus it would seem 
that only a highly contrived and flnc-tuncd cosmology could exhibit the non-isotropizing 
property and consistency with the data. However, during the last decade, as mentioned 
in the introduction, evidence has accumulated against both of the presumptions above. 
Therefore, it is interesting to look for fixed points, especially scaling solutions, in the more 
general set-up allowing for direction-dependent expansion rates. It becomes then possible 
to see which kind of possibilities in general exist for anisotropic expansion histories, and 
whether they could occur naturally (i.e. without fine-tuning). We therefore perform a 
dynamical system analysis. Then, we will also consider the specific case of a generalized A 
term. 
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2.1. Parametrization 

Consider a general fluid flow in a curved space-time. Assume the 4-velocity tangent to 
the fluid flow hues is u". This is the central object in the formulations below. Everything 
that follows will be defined in relation to this average flow. An useful tensor is then the 
projection tensor orthogonal to m", which is given by 

hab = gab + UaUb- (1) 

The time derivative of any tensor T"^ -^^ is understood as the derivative along the vector 
u"-, as 

fab...^^ = T^'-^,,,^. (2) 

Note that in general this differs from the derivative with respect to cosmic time, and agrees 
only for scalar fields. We also define ^^^V, the covariant derivative operator obtained by 
projecting the four-dimensional one as 

VeT'"''"'^^ = N^g^ih^y ■ ■ - hf hf . . .T"''' "c'd'...;e'- (3) 

The first covariant derivative of the u"' is decomposed as follows: 

Ua-b = + (^ab + ^dhab " Oa^^b- (4) 

One calls Oa = iia the acceleration, 9 = u^.^^ the expansion scalar, ujab — <^[a6] the vorticity 
tensor and (Tab — cr{ab) the shear tensor. Both the shear tensor and vorticity tensor are 
traceless and orthogonal to u^. It is useful to define a representative length scale i by 

This is a covariant generalization of the usual e-folding. Furthermore, the stress-energy 
tensor of a general fiuid in curved spacetime is decomposed as follows: 

Tab = PUaUb + phab + "^QiaUb) + T^ab, (6) 

where p is the pressure, p is the energy density, qa is the energy fiux and -Kab the fiuid 
viscosity. 

In the following we consider the choice u"- — 5^, so that the fiuid fiow is given by the 
time-like unit vector. We also restrict ourselves to homogeneous cosmologies, and thus set 
any spatial gradient to zero. Furthermore, we assume that the vorticity vanishes. These 
assumptions, summarized as 

u'^^l (3)V„/ = 0, a;„6 = (7) 

characterize the Bianchi I cosmology. If we would further set Uab — T^ab = 0, we would 
recover the Priedmann-Lemaitre-Robertson- Walker (FLRW) universe. 

In the Bianchi I spacetime, one needs three parameters to characterize the pressure of 
a fluid in all directions. In addition to the isotropic part which one gets as the trace part 
of Eq.(6), 

P = Ih'^'Ta,, (8) 
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pia, = hlh% - -h'"^ T^. (9) 



two more parameters are needed to fully determine the degrees of freedom in the traceless 
and symmetric tensor 

1 

''ia"'b) - 3' 

In Bianchi spacetime one may choose diagonal coordinate system, as we'll discuss later. 
It is occasionally claimed that one may consider the axisymmetric case without loss 
generality. It might then seem that one always choose the anisotropy axis parallel to 
one of the coordinate axes and thus, by considering a coordinate where anisotropy along 
one direction vanishes and due to the tracelessness requirement must in the remaining 
directions have an opposite sign and equal magnitude, describe the anisotropic pressure 
with full generality. However, rotation of the anisotropy axis would not preserve the 
diagonality of the system. The 3x3 traceless symmetric tensor in (9) has 5 independent 
entries, but 3 of them may be fixed by choosing the spatial coordinate system. Thus 
the anisotropic pressure is characterized by two degrees of freedom. This resembles the 
anisotropic property of gravitational waves, which, at the level of linear perturbations 
are conventionally parametrized by the two modes of polarization with respect to the 
wavevector of perturbations. Since in the homogeneous background there is no obvious 
definition for such a wavevector, we instead use the following contractions for a fully 
covariant characterization of our anisotropic degrees of freedom: 

L»=lw<^^ T=l7r^>^ (10) 

Higher contractions would not be independent of these, but instead Tr^Tr^Tr^vr^ = 2D^, 
'^b'^d^a^c^e — 1^-^) -'■^ following wc employ a combination of the quantities in 
(10) as 



where 



\ D^F 5 / , 



F^^D^- 12T^ + 2V6V6r4 - D^T^. (12) 



Prom these covariant expressions the use of this parametrization is far from obvious. Let us 
therefore make an excursion to the preferred coordinate system picked up by the symmetries 
of this model. In the diagonal system, where tt = (tt^,, tt^, tt^), we have now 

S ^ - ^{7rx-7ry)/px, (13) 

7 = - ^(27r^ + 7^y)/px■ 
The use of these variables is mainly motivated by their straightforward interpretation as 
generalized equations of state (see the Appendix) and concordance with previous notations. 
The parameters 5 and 7 may be interpreted simply as differences of the pressure along the 
X and the pressure along the y and z axes, respectively. The axisymmetric cases are now 
seen to correspond to the three possible cases that either 5 = 0, 7 = 0, or 5 = 7. One 
might get rid of the ambiguity by relabeling the axes in such a way if the pressure along 
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any two axes is equal, the third one is called x (a general rotation is not of course allowed 
by the symmetries) . In the similar way, we will parametrize the shear stress of the metric 
as: 

3 

R=-Q{(^x-(Jy). (14) 

Again, this is just a change of variables which can be convenient for some purposes. 
The dimensionless variables R and S might be interpreted as the fractional difference of 
expansion rates between the y and x (for Ft) and between the z and x (for S) . The notations 
(13) and (14) have been introduced in Ref. [49]. Again, the straightforward interpretation 
of the shears in the metric sector, R and S, is confined to the preferred coordinate system, 
but these quantities can be defined fully covariantly in a way completely analogous to 
coefficients of the shear in matter sector, 5 and 7 in Eq.(lO). 



2.2. The matter content 

As the contents of the Universe we consider a two-fluid system, where only one of the fluids 
is responsible for possible nonstandard (anisotropic) properties. We have then a perfect 
fluid with an energy- momentum tensor (we denote it with index m) and an imperfect 
energy-momentum tensor (this wc denote with index X) 

Toh = Pm{UaUb + WmKb) + Px{UaUb + Wxhab) + T^ab, (15) 

where Wm is the equation of state parameter, Wm = 1/3 for radiation and Wm = for 
dust (and something in between for their mixture). We also consider an imperfect fluid, 
which we allow to have the most general energy-momentum tensor compatible with the 
assumptions (7). The shear stress Hab is due to this component solely, and therefore we 
drop the X from this quantity. 

The continuity equations are then given by the divergence of the energy-momentum 
tensors. We let the two components also interact, and thus allow nonzero divergence for 
the individual components. 

Pm + (1 + Wm)Opm = QPm, (16) 
PX + (1 + Wx)9px + (TabT^"^ = -QPm- (17) 

where Q describes the coupling. If there are no interactions, one finds that the matter 
density scales as Pm ~ (i^)^^^'^"^. For most of the study, we will simply neglect the 
coupling between the X and the isotropic component. However, we keep it for the sake of 
generality, since wc will later be interested in flnding cosmological scaling solutions. 

The system is of course not closed until we have determined the properties of the 
imperfect stresses TTf^u- This is equivalent to specifying the anisotropy parameters 5 and 
7 in Eq.(13). In the case of viscous fluids [62], these stresses are related to the expansion 
factors. A covariant form for the viscosity generated in the fluid flow is 

TTab = (ua;ch\ + Ub-ch\ - ^H^.Jlab] + CW.Jlab- (18) 
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Model 


5 


7 


Wx 


Isotropic fluid 








Wx 


String 





1 


1/3 


Domain wall 


1 


1 


2/3 


Magnetic field 





-2 


1/3 


Viscous fiuid 


9p 


9p 




A, case I 


\e{2S - R) 


^0{S - 2R) 


-1 


A, case II 


5 (free) 


7 (free) 


-1 - 1(3 - 2i? + S)5 - 1(3 + i? - 2,5)7 



Table 1. Models with anisotropic equations of state. 



Now the conservation equations T'^.^ — reduce to the Navier-Stokes equations in the non- 
relativistic limit. Here q is the shear viscosity coefficient, and C represents bulk viscosity. 
Heuristically, they describe the resistance of fiuid fiow to external stresses. Such result in 
acceleration (derivatives of the velocity four- vector) which are then compensated by effect 
of TTab- The bulk viscous part reacts to volume- changing stresses as it is proportional to the 
expansion scalar; in cosmology it gives an extra Hubble friction. The shear viscous part 
corresponds to the the symmetrized and traceless part of the acceleration, thus reacting 
to an anisotropic stress. Ultimately the coefficients C and q should be determined from 
the microscopic kinetic theory. Calculating the tensor (18) explicitly with the help of the 
covariant definitions of section (5.1), one can write the equations of state as 

w^wx-[3C + ^^iR + S),]-, 5^^, 7=^. (19) 

One notes even with nonzero shear viscous coefficient q, one does not generate shear in 
the metric. In fact, with a positive coefficient q (which is required by second law of 
thermodynamics), the viscous property of the fiuid in fact tends to decrease the amount 
shear J. A bulk viscosity [63-66] could also have interesting consequences in this framework; 
however we do not focus further on the Navier-Stokcs type of stress in the present study. 

There are also several examples of imperfect cosmological sources with constant 
equations of state. A magnetic field along the 2;-direction [67] could be modeled by 
Wx = 1/3, (5 = 0, 7 = —2. A string with a constant tension along the 2;-direction [68] could 
be described by wx = 1/3, 5 = 0, 7=1. An anisotropic cosmological constant is described 
by Wx = — 1- We mention these possibilities in Table 1. We will return to these cases later. 
Note that our description with constant skewness parameters unifies and generalizes all of 
the mentioned models. Vector fields will be considered in much more detail in elsewhere 
[22]; we have found that such models can also feature constant equations of state S and 7 
in special cases, though they in general are dynamical, depending on the potentials and 
couplings [22]. We have gathered some other examples of imperfect matter in Table 1. 



I This may be seen by plugging the S and 7 in the evolution equations (25,26). 
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2.3. Evolution equations 

The generalized Priedmann equation can be written as 

^2 = 3(p^ + px + <j'). (20) 
We also define the dimensionless density fractions 

n„^ 247rG^, U^-P^. (21) 

Pm + PX 

We introduce Jl^ to quantify the amount of matter in analogy to the standard cosmology 

as the percentage of matter contribution to the square of the average expansion rate. 
Note that there is an ambiguity when referring to matter densities since, the generalized 
Friedmann Eq. (20) implying 

Q„ = (l-C/)[1- 3(77^2], 

Vtm is different from the matter density divided by the total density (which is 1 — f/). 

Our dynamical variables are then the density fraction f/, and the two shear anisotropics 
R and S. As the time variable we define x = \og£. In the following a prime denotes 
derivative with respect to that. The evolution equations for these are then 

U'^U{U-1) [n'^'aab + 3{wx- wj] - UQ, (22) 
R'^-[3 + jjR + 96{—-l)U, (23) 

'^' = -(3 + ^)^ + 97(^-1)^/. (24) 

One notes that only with non-vanishing skewness parameters, shear could be generated 
also in an initially isotropic situation. Also, ii 5 — ^ — 0, even initially nonzero R or 
S would typically decay. Although the coupling term Q does not comphcate this system 
much, but appears only in the evolution equation for U, its presence of course can change 
the dynamics completely. As mentioned earlier, we consider mainly Q = case in this 
paper, but keep the coupling in the equations for the sake of generality The two last 
equations can be written explicitly in terms of solely R, S, and U as 

^'^'^ ~ "^"^^ +Wm-1]- 67C/} (25) 
i?' = ^ (1 - ^) {R [U {wx - wm) +wm-l]- 6SU} (26) 



In the following we will refer to an effective equation of state, which we define as 

2 9' 
1,1 

if (T^ = 0, the universe expands as if dominated by a fluid with Wm = Weff] ff cr^ 7^ 0, the 
average of the three effective equations of state in each direction is Wg//. 



^eff = ~- - 1, (27) 
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3. The Background as a Dynamical System 

3.1. Fixed Points in the Axisymmetric Case 

In the Friedmann-Lemaitre universe it has been proven difficult to address the coincidence 
problem by finding a model entering from a matter-dominated scaling solution to an 
accelerating scaling solution. Allowing the presence of three expansion rates perhaps opens 
up the possibility of finding more such scaling solutions which might eventually help to 
understand the coincidence problem. It is also possible to model infiation within this 
context (an interesting case then describes inflation as an era of primordial isotropization 
to the solution R = S = from a more general initial state). Therefore we consider the 
asymptotic behaviours of the universe. In the simplest anisotropic case we may assume: 
1) axisymmctry, 5* = 7 = 0, 2) constant equations of state and skcwncss parameter 
Wm = w = (5 = 0, and to begin with, also 3) no coupling between the components Q = 0. 
To end the analysis, we briefly review the consequences of relaxing the assumption 3). The 
results are summarized in the Table (2), and also discussed in some detail below. Please 
note that now R can be considered as a shorthand notation for the expansion-normalized 
invariant shear 

R = 3v^^, 

u 

and similarly 5 corresponds to the density-normalized invariant shear 



Our system has 5 fixed points. With the assumptions l)-3), there are two fixed points 
corresponding to matter domination. One is the 

• FLRW solution 

i? = 0, [/ = 0. (28) 

We assume Wm < 1- Then the FLRW solution is a stable node when w + S > Wjn, 
otherwise a saddle point§. 

• An anisotropically expanding empty universe 

i? = ±3, [/ = 0, (29) 

is also a fixed point. This is a saddle point in the R = +3 case when 6 < w — Wm, 
and in the R = —3 case when 6 > {wm — w)/3 and otherwise unstable. This is the so 
called Kasner solution, where the metric may be described by 

d^Kasner = 'dt' + t^^^dx^ + e^^dy^ + ^^^dz^ (30) 

where the Kasner exponents satisfy 

Ka + + K,^ Kl + Kl + Kl (31) 

§ In this section w corresponds to the equation of state in the x- and y directions. Thus wx = w + 5,].i 
wx stands for the average of the principal pressures. 
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In the present case, due to the axisymmetry, the exponents Ka = are of course 
equal, but we wrote the general form for later convenience. This solution is relevant 
near the singularity, t — >^ 0, where the matter may be neglected, Pm/ H"^ ~ p/H^ ~ 0. 
This solution may be realized in two ways||. Firstly, one may have expansion in two 
and contraction in one direction, so that Ka — —1/3 and K^, — Kc — 2/3 (or let 
the two directions change roles a <-> 6 to get the solution with R — —3). Secondly, 
one may have two static directions and one expanding direction, corresponding to 
Ka — —1/3 and Kf, = Kc = (or again let the two directions change roles a b to 
get the solution with R = —3). 

• An anisotropically expanding universe dominated by dark energy, 

= ±3, [/ = 1, (32) 

is also a fixed point of our system. The pancake solution with i? = +3 is a saddle 
point ifw>l + SoYw< Wm + S, and the cigar solution with i? = — 3 is a saddle 
point if w + 35 > 1 or w + 3(5 < Wm- Otherwise they are unstable. On the average, 
the equation of state is for a stiff fluid, Weff = 1. 

• The anisotropic fixed point 

R^ C/ = l. (33) 

represents a balanced expansion in different directions. Assuming here for simplicity 
that Wm — and restricting to the case w < Wm, this can be a stable node if either 
5 + w > 1 or 

1 + 66+ ^l + 16w{w- 1) > 2w A 2w + ^/l + 16w{w - 1) > 1 + 66. 
We have a saddle point if either 

36 + w <lA2w + yi + 16w{w - 1) <l + 66 

or 



l + 6>wAl + 66+ + 16w{w - 1) < 2w. 



The scaling solution 



^ ^ 3{w + 6 - Wm) ^ {w + 6 - Wm) {Wm - 1) ^g^^ 

26 ' 3(5^ — 26 {w — Wm) — (w — Wm)^ 



is finally the fifth fixed point of our system. Considering again only cases when 
w < Wm — 0, this is a saddle point if either 

5 > A 2w + + 16w{w - 1) >l + 66 



or 



6 <QM + 66 + yi + 16w{w - 1) > 2w. 



General stability conditions for this point are more complicated. Instead of trying to 
find feasible analytic expressions for them, we resort to an empirical study. 

II In terms of the expansion rates introduced in the appendix, tlic mentioned two possibilities correspond 
to d/a = — l/(3t), and b/b = c/c= 2/{3t), and a/a = 1/t, and b/b = c/c = 0. 
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To study when the solutions are relevant for the dark energy problem, we use a dust- 
dominated FLRW universe, Wm = R = S = 0, a.s the initial condition. Then the outcome is 
uniquely determined by dark energy properties, w and 5. The numerical results, shown in 
Figs. (1,2) and summarized in Fig. 3, agree completely with our analytical considerations. 




Figure 1. Phase portraits in the cases S — —0.25 (right) and 6 = 0.25 (left) when 
w = —1. The solution (28) is at i? = [/ = 0, and the fixed points (29) and (32) arc at the 
four corners of the portraits. In both cases the fixed point (33) attracts trajectories from 
everywhere. 




Figure 2. Phase portraits in the cases S = —2 (right) and ^ = 2 (left) when w = —1. The 

solution (28) is at i? = J7 = 0, and the fixed points (29) and (32) are at the four corners 
of the portraits. In right panel all trajectories lead to the scaling solution all (34). The 
left panel shows isotropizations in the example with w + S > Wm- 



Let us finally consider relaxing the assumption 3) to include a coupling Q. We assume 
it is constant. This may be motivated by the analogy with scalar fields (exponential 
interaction term in the Lagrangian results in constant Q for scaling solutions) , and can also 
be realized with vector fields having matter couplings, as we show in another publication 
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Figure 3. The asymptotic state of tlie universe starting from an Einstein-deSitter stage. 
Tire future fate depends on the dark energy properties w and 5 and is classified into three 
possibihties. Ii5+w < 0, the isotropically expanding dust domination (solution (28) where 
U = 0) continues forever. Otherwise, the universe will end up expanding anisotropically 
and either dominated by dark energy (solution (33) where [/ = 1) or exhibiting a scaling 
property (solution (34) where < J7 < 1). Fig. (1) showed phase portraits of two U — 1 
models, the right panel of Fig. (2) showed an example of the < J7 < 1 case, and 
finally the left panel of Fig. (2) showed isotropizations in the U = (here we refer to the 
asymptotic value of [/). 



with explicit examples [22]. Then the dark energy dominated anisotropic solution Eq.(32) 
would turn into an anisotropic scaling solution 

R = ±3, U = l-- (35) 

3(w„ - w) + (3 =F 6)5 

and where the amount of matter is proportional to the strength of the coupling Q. 

The scaling solution (34) generalizes in the presence of a coupling to 

15(^2 + (_3 + g + sw){w -Wm) + 6{Q + Qw- 3(1 + «;„)) ± ,/? 



R 



4:6{-l + 6 + w) 



jj _ - 6{3 + Q + Qw - 9wm) - {3 + Q + 3w - Qwm)iw - Wm) ± v^ (^3g^ 

6(35 + w - Wm)i5 - W + Wm) 



where 



r = (—95^ + 5(3 + Q + 6w — 9wm) + {3 + Q + 3w — 6wm){w — Wm)) — 

12(35 + Q + 3W — 3Wm){36 + W — Wm){Wm — 1)(5 — W + Wm)- (37) 

Thus the coupling modifies both the relative amount of anisotropy and the relative amount 
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(32) 
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±3 
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See Eq.(35) 


DE 2 


(33) 


1 4<5^ 

5+w-l 


w+6-1 


1 




Scaling 


(34) 


W \ "I ) 


3(«) + 5-«)m) 


{w+5-Wm.){Wm-i) 


See Eq.(36) 




25 


35^ — 2(5(ui-UIm) — (f — Ulm)^ 



Table 2. Fixed points in the axisymmetric case, with constant equations of state w, 
Wm and S. The effective equation of state We// is defined in Eq.(27). The last column 
indicates whether the solution can be modified by a coupling. The fixed point corresponds 
to an anisotropic expansion unless R = 0. 



of matter. An isotropic scaling solution is now possible with a nonzero coupling and 5 — 0. 
We list a summary of the fixed points found in this subsection in the Table (2). 

The rest of the fixed points, Eqs. (28), (29) and (34) would retain their form in the 
presence of a coupling term. 

3.2. An anisotropic Cosmological Constant 

We will show that there does not exist an anisotropic generahzation of the cosmological 
constant with constant equations of state. Then we consider two different classes of time- 
dependent cosmological terms with constant energy density. 

The continuity equation (93) for a component with general w, 5 and 7 reads 

PA + eil + WA)pA + 7r„6(T«^ = (38) 



PA + 



;i + wa) + Ij{R- 2S) + y {-2R + S) 



PA = 0. (39) 



We have added the subscript to emphasize we are discussing a cosmological constant term. 
In the second line we wrote open the component contributions to the invariant shear. As 
we then want /?a = 0, the square brackets term should vanish identically. If we have 
time-independent equations of state, we cannot consistently model an anisotropic A term. 
Obviously, the constant terms in the square brackets should then vanish, and we should 
set wx = —1- Since there also appears R and S, which in general depend on time, the 
linear combination 

X = ^7r„,a"^ = {R- 2S)j + {-2R + S)5 (40) 
u 

should however be zero. To see how this is impossible, one derives an evolution equation 
for X: 



e+^-\x + 3 [7(5 - 27) + 5(7 - 26)] (^ - l)eU. (42) 
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The second term will always enforce X to evolve unless 5 = 7 = 0. Note that the second 
term does not vanish since Eqs. (25,26) set the shears evolving. 

As a summary, we note that in terms of our covariant notation, the two conditions 
mentioned above correspond to the two statements 

• wa = 1, i.e. the isotropic part of the pressure should be equal to minus the energy 
density, pa = -pA- 

• X = 0, i.e. the anisotropic part of the pressure should be orthogonal to the shear 
expansion, Tr^^Uab = 0. 

If both of these conditions are satisfied, with constant 5 and A, these must be zero. We will 
therefore consider dynamical equations of state in the following subsection 3.2.1. However, 
it turns out that to find interesting cosmology one may have to violate one of the conditions 
above too. We will consider this in the later subsection 3.2.2. 



3.2.1. Case I: Constant Isotropic Pressure. One may satisfy both of the two requirements 
stated above with time-varying skewness parameters. We introduce a parameter h and 
write 

5 = ^e{2S-R), (43) 

^^^e{S-2R). (44) 
Then the anisotropic stress becomes 

7r% = diag [0, -S + R, S, -R] hH, (45) 
which indeed is, by construction, orthogonal to the metric shear 

a\ = diag [0. R + S, -2R + S,R- 2S] H/3. (46) 

Since now also Wx = — 1, both the energy density and the isotropic part of the pressure are 
constant. This prescription seems to be equivalent to the model introduced by Rodrigues 
[53], who also shows it is possible to associate such properties of the cosmological term 
with an infrared non-commutative property of the spacetime. The anisotropy arises from 
the deformation of the Poisson structure between the canonical momenta in such a way 
that the dynamics remain invariant under rescalings of the scale factors [53]. The non- 
commutativity coefficients are then proportional to the parameter h (which in our case has 
dimension of 1/^) . 

This would be difficult to apply for description of anisotropics in the present universe. 
The reason can be seen from skewness parameters (43) and (44). Since they are 
proportional to the shear, the model does not support spontaneous generation of anisotropy. 
If R and S vanish in the early universe, they will stay zero even when the A term becomes 
dominant. Observable deviation from ACDM could occur if the system turned out to be 
unstable: the A term could in principle react to even perturbatively small anisotropies by 
amplifying them into direction-dependent background expansion. Let us briefly study this 
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possibility. After algebraic manipulation one may write the evolution equations as follows 
U'^ -W{U -l){l + wm), (47) 
E!^^{l- [-mSU -R{l + U- 2h0U + {U- l)wm)] , 

S' ^^(^1- [4h0RU -S{l + U + 2h9U + {U- l)wm)] , 
where the Hubble factor can be expressed as 

= ^ , (48) 

[t^(l-^)]' 

with the rescaled constant h which is now dimensionless, defined as h = hy/2A7rGpA- The 
equation for U' is of the same form as in ACDM cosmology. The solution is thus 

= l + (_^_l)e-3(i+«.^)x' 

where Uq is the relative amount of vacuum density today. It is also obvious that R — S — 
is a solution with any value of w. In fact the system has several fixed points. 

• The Kasner solutions 



R = ^±^9-ls^. (50) 

These correspond to the metric (30). Since we are not restricted here to a specific 
axisymmetric case, the Kasner exponents have only to satisfy the relations (31). The 
solutions with different — 3 < S* < 2 correspond to these different exponents. Since for 
the Kasner solution matter is negligible, the value of U is irrelevant (though formally 
should be set either to U = or U = 1). One quickly notices that this class of solutions 
is unstable: the anisotropy decays and matter becomes eventually dominant. 

• The A dominated isotropic solution 

R^S^O, U = l. (51) 

If we then consider small perturbations U — 1 + u, R — + r, S — + s about the A 
dominated solution, we find that 

u ~ x-^^^+'""^\ while r, s ~ x'^^'^'^. 

Hence, the anisotropy decays, proportional to the inverse volume element. If 
anisotropy is inserted as an initial condition, the non-zero skewness parameters (43) 
and (44) only cause an oscillation about the usual diluting behaviour. 

• The matter dominated isotropic solution 

R^S^O, [/ = 0. (52) 
One may also check that for the matter dominated solution, the perturbations go as 
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The negativity of the first exponent only signals that the fractional matter density is 
diminishing. The two last ones tell us that during matter domination the anisotropy 
decays too, though slower than during the A domination, and to first order without 
oscillations. 

For a cosmological evolution the relevant solutions are then the two last ones with 
[/ = and the other with U — 1, and both having R — S — These will come 
out according to Eq.(49). The anisotropy is neghgible by definition when the A term 
can be neglected. Shear is not generated during the evolution of the universe, since 
anisotropic perturbations about the usual solution decay. We have checked numerically 
that this behaviour is general and is not restricted to the asymptotic regimes of cosmological 
evolution. 



3.2.2. Case II: Constant Anisotropic Pressure. Let us now take a different approach, and 
keep 5 and 7 as free parameters describing the anisotropy of the A term. The Eq.(38) 
should still reduce to the identity p\ — 0. There is a simple prescription which achieves 
this by forcing the isotropic pressure to compensate for the anisotropic part: 

WA = -1- ^X. (53) 

The equations of state associated to the A-tcrm can now be also different from —1. Also 
in such cases the vacuum energy density stays constant, but then the pressure becomes 
direction dependent, forcing eventually the universe to expand anisotropically when the 
cosmological term is significantly large. The evolution equations for the background 
variables then reduce to 

U'^ -3U{U-l){l + Wm), (54) 
S' ^^(9-R^ + RS-S'^)i^S -U (^1 + ^X + +Wm-1 - 67C/| 

R' ^^{9-R^ + RS-S^)i^R -U (^1 + ^X + +Wm-1 - GSU^ 

The equation for U' in (54) is of the same form as in ACDM cosmology, and is given by 
Eq.(49). We have not found an analytic solution for the shears it! and S. 

For simplicity, let us consider the axisymmetric case that 5" = 7 = 0. We then find 
several fixed points. 

• The FLRW, 

i? = 0, U = 0, (55) 

corresponding now to Eq.(28). It is easy to sec that this point is always unstable. The 
matter domination will inevitably be followed by the cosmological constant dominated 
era, which is anisotropic whenever 5^0. 

• An anisotropic matter domination, 

= ±3, U = 0. (56) 
For dust matter this is always an unstable point. 
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• An anisotropic A domination I, 

= ±3, [/ = 1. (57) 

For dust matter, this is a saddle point. 

• An anisotropic A domination II, 

R = ^(^± \/9 + 36(52 j ^ U = 1. (58) 

One may show that when Wm = 0, this is always an attractor. More specifically, 
introducing small perturbations about this solution, and constructing the matrix for 
the derivatives of these perturbations from the system (54), the eigenvalues g„, qr of 
the matrix then are 

?« = -3, qr = -^VlTW(±l + VlTWy (59) 

Since these are negative regardless of the value of 5, we know that the perturbations 
tend to decay. 

A scenario then appears naturally, featuring a transition from the usual isotropic fixed 
point (56) to the anisotropic vacuum dominated point (58). This point coincides with the 
solution (57) when 5 is very large. One notices that the asymptotic value of the shear R 
goes to R ±3 when 6 ±oo. We have numerically verified that the universe always 
evolves along this track when begun from an isotropic initial stage. This then provides the 
possibility to generalize the standard ACDM cosmology in such a way that the universe 
features anisotropics at late times when the A-term begins to dominate. 

If the matter content is dust, the equation of state for the total fluid in this model is 
given by 

WTOT = -(1 + j)U, (60) 

and the total sound speed squared is 

2 1 \ _ 5(i?2-9)C/[185C/ + i?(3 + (3-5(3 + 2i?))C/)] 

CsTOT - 9^ l^x^T ) - 81(C/-1) '^^^^ 

assuming that S is constant. One notes that this could be positive or negative. It vanishes 
when S does: the usual ACDM universe is "silent" in the sense that the CDM is pressureless, 
no perturbations propagate in the (isotropic) A medium. The anisotropy quantity ^ is now 

"^HR 

It was expected that this is proportional to the shear in the expansion rate, HR. 

In the Figure 4 we show the asymptotic state of the universe as a function of the 
skewness parameter 5. In this numerically produced plot we also consider some cases 
where 7 does not equal zero, (though for simplicity we restricted to an axisymmetric 
case 5" = 7 = in the analytic considerations). We return later to briefly consider the 
observational implications of this anisotropic generalization of A. 
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Figure 4. The asymptotic state of tlie universe starting from an Einstein-deSitter stage 
for tlie anisotropic generalization of cosmological constant having constant 5 and 7 (case 
II, section 3.2.2). The R (blue, solid lines) and S (red, dashed lines) are given as function 
of 5. The thick lines depict the cases 7 = 1,2,3,4,5, and the thin lines depict the case 
that 7 = 0,-1, —2, —3, —4, —5. The analytic solution (58) is recovered in the case 7 = 0. 
These curves show how this generalizes non-axisymmetric cases. 



4. Observational Bounds from an Anisotropic Background Expansion 

The part of this study focusing on homogeneous cosmology is ended by a comparison 
of the background dynamics of these models to the cosmological observations. In 
Section (2) we have made a phenomenological description of an anisotropic component 
in terms of the skewness parameters 5 and 7, Eq.(13). The question of its observational 
consequences is taken under investigation now. We then focus on dark energy and an 
anisotropic background expansion. Observable effects which might have been imprinted 
during the early inflationary era, would appear to be more straightforwardly described as 
perturbations and thus not changing the FLRW predictions at zero order [69,70]. As we 
emphasise later, a proper treatment should take into account perturbations, in particular 
(the possibly anisotropic) initial conditions set at inflation [71]. However, as an attempt 
at a model-independent first constraint on the model, it is useful to estimate the bounds 
ensuing from the background expansion. 

4-1. CMB Anisotropies 

A possible anisotropy will leave its imprint in the CMB. To calculate the CMB spectrum 
in a background with a Bianchi I metric Eq.(122), we begin from the geodesic equation for 
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photons, 

where the only nonvanishing connection coefficients are (summations are suppressed for 
Latin indices in this section) 

r°=a,a„ r^, = no = ^. (64) 

We are interested in the redshift of photons coming from different directions. The four 
vector — dx^ j dX is subject to the null condition, 

«% = = -(^)%^a,^(pf , (65) 

where we call = P\ Consider two photons, the ffist emitted at t = and the other a 
tiny bit r later at t = te + t. The difference of the two null conditions gives 

= E«^«^ (^f ^(A) + 0{t'). (66) 

i 

The definition of redshift is 

l + .(Ae)^^, (67) 

where r(Ar) is the time difference of the received signals. Differentiating Eq.(67) and using 
Eq.(66) yields 

dlog{l + z) 1 . . , . 

i 

The P* can be obtained from the spatial components of the geodesic equation (63), 

^ + 2^P'k° = 0. (69) 
dX tti 

We find P* ~ Oj"^. Since one can always rescale the scale factors, let us fix ai(to) = 1 
where to indicates today. We then write our initial P* as P*(to) = Pi- Since one 
can reparameterize A without affecting the results, we choose now such a normalization 
that X^jP? — 1- t^^s think about the unit vector p in terms of the angles 

{Px,Py,Pz) — (sin ^ COS 0, sin ^ sin 0, COS ^). Thus we have managed to parametrize the 
photon paths according to the angle they hit the observer today. Inserting this solution 
for P* and from Eq.(65) into the evolution equation for the redshift (68), one finds that 
it can be solved to yield 




1 + ^(P)=I>.^I • (70) 
Since the sum is equal to 



a • 2 J. 1 ^ 2 



sin'' 9 I cos + 7^ sin"' + — cos^ 9 



(71) 
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it is easy to see that one can rewrite the result as 

l + z(p) = -Jl+p2e2+p2e2 (72) 

(X V 



in terms of the eccentricities 



The temperature field is determined by this redshift by the relation 
and the spatial average is 

47rf = J dnpT{f)). 

The last scattering temperature does not depend on the direction, but the photons 
coming from different directions will be redshifted different amounts. The anisotropy field 
is then 

^^1-^. (75, 

The coefficients in the spherical expansion of this anisotropy field are called a^^, and due 
to orthogonality of spherical harmonics Y^^, are given by 

aem = I dnJ^Yl^. (76) 
The multipole spectrum can be described by 



1 £(£ + 1) ^ 

2^ \aimV. (77) 



^ 27r (2£ + 1) 

One may expand the redshifts (72) in the eccentricities (73). Then one notes that the aim 
will be real (since there is only even dependence on the polar angle in the anisotropy field 
and the imaginary parts of the integrate to zero) , and that for all odd I the a^m will 
vanish (since only even powers of the azimuthal appear in the expansion). To first order 



6 

and in addition to the monopole, there is only the quadrupole 



1 pK I 2 2\ 

which implies 



021 = 02-1 = 0, a22 = 02-2 = - \\ -^^y 




2 



Q2^^^e%-re\-elel (78) 

The observed value of this is Qiiphs) ^ 5.8 10~^, while the standard concordance model 
predicts Q2{iso) 1.2985 10~^. It has been suggested in previous works also that this 
discrepancy could be explained by an cllipsoidality of the universe [51,52]. This would 
require that the anisotropy of the background is suitably oriented with respect to the 



CONTENTS 



23 



intrinsic quadrupole and cancels its power to sufficient amount. For any orientation then, 
we should have Q2 < 2.7209 10^^ to be consistent with observations taking into account 
the cosmic variance. The constrains this implies on the dark energy equation of state 
parameters are rather tight. 

One should, however, keep in mind that we have only studied the Bianchi I type 
background, which leaves a very simple pattern on the CMB, and assumed that to be 
somehow superposed with the anisotropics of the standard perturbed FLRW universe. 
This is of course a gross description. In a more realistic treatment, one should study 
the perturbations in matter in the anisotropic background. This is not a straightforward 
problem and will be tackled in section 5. We now focus only on the background anisotropy 
and assume that at least the order of magnitude of the limits one derives from it reflect 
correctly the properties of the universe. 

In more general models, e.g. with time-varying 5 and 7, one could allow more 
anisotropy. It is in principle possible for arbitrarily anisotropic expansion to escape 
detection from CMB as long as the expansion rates evolve in such a way that Cz — Cy — 
In other words, the quadrupole vanishes, if each scale factor has expanded - no matter how 
anisotropically - the same amount since the last scattering ^. 

Generating these effects at low redshift has an advantage that it relaxes constraints 
which would otherwise come from the CMB polarization [69] and could be strong for a 
given temperature anisotropy in isotropizing models because of the signiflcant polarization 
anisotropy at last scattering. However, anisotropic dark energy could evade this since the 
optical depth to 2; ~ 1 is very small. More specifically, polarization of the CMB photons 
is generated by the Thomson scattering process on the free electrons. When the electron 
distribution has quadropole anisotropy, the scattered radiation will have net polarization. 
Thus it is clear that an anisotropy of the expansion will in principle always have imprints 
on the polarization pattern. However, most of the photons we see in the CMB were last 
scattered at z ~ 1100. Any model which predicts large anisotropics at such redshifts could 
thus easily be ruled out. At smaller redshifts, the effects of anisotropy on polarization 
are only secondary, since the photons do not scatter, but their energy distribution may be 
slightly modified. At the reionization epoch one expects again some photons to scatter, 
but the impact this has on CMB is restricted to the largest scales and suppressed by 
several orders magnitude. Thus it seems easier to explain the observed anomalies in the 
temperature spectrum without running into problems with unobserved anomalies in the 
polarization spectrum by anisotropies occurring in the later rather than in the earlier 
universe. 

4.-2. SNIa Luminosities 

Anisotropy in the acceleration rates would be in principle seen as rotationally non-invariant 
luminosity distance - redshift relationships for the SNIa. The main CMB constraints on the 

^ We emphasize that this is only true for the CMB pattern due to the background anisotropy. In the 
next Section of the article we consider inhomogeneous perturbations, which in general have anisotropic 
signatures also in the case that they vanish at the background level. 
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present model ensue from the quadrupole moment subject to effects of cosmic variance, and 
at least in specific models, the constraints can be loose or nonexistent. It is then interesting 
to see if one could limit the anisotropy with the present data on the nearby supernovae. 
More accurate SNIa data will become available in the near future, notably from the SNAP 
experiment"*", but we find that useful constraints can already be derived. The luminosity- 
redshift relationship of the SNIa can be used to probe the possible anisotropies in the 
expansion history. This is a complementary probe to the CMB quadrupole, since the SNIa 
objects arc observed at the z < 2 region, whereas CMB comes from much further away 
at z ~ 1000. The luminosity distance as a probe of anisotropies has been considered in 



The spatial geodesic eq.(69) tells us that the direction p of a photon that is coming 
towards us is constant. This reflects the orthogonality property of the Bianchi I model, in 
contrast to tilted models where the geodesies would involve rotation effects. For the null 
geodesies we have then 



To test this prediction with the data, we apply the formula (72) for each observed rcdshift 
of a supernova and match its distance modulus }A inferred from the observation to the 
one computed from Eq. (80) by 



We have to then take into account also the angular coordinates of each individual 
supernovae in the sky which fix p for each object. In our analysis we use the GOLD data 
set [74], which consists of five subsets of data*. We marginahze over the directions in the 
sky by integrating over the likelihood-weighted over three Euler angles and normalizing. 
Similarly we also marginalize over the present value of the Hubble constant. Taking these 
cosmologically irrelevant parameters into account makes the computations much heavier 
but is necessary. The results are summarized in Fig. 5. The best-fit anisotropic models 
are only slightly preferred over the ACDM, the difference being A^^ ^ 1. On the other 
hand, the SNIa data allows skewness in dark energy to an interesting degree. 

In Fig. 6 we illustrate the SNIa constraints in the axisymmetric cases. Larger skewness 
\5\ would typically be compatible with the SNIa data for equations of state w < —1 and 
large matter densities. This means that even if the CMB formed isotropically at early 
time, it could be distorted by the acceleration of the later universe in such a way that it 
appears to us anomalous at the largest scales. The future SNIa data, with considerably 

+ see http://snap.lbl.gov/. 

* We find the angular coordinates of each of the 182 GOLD supernovae partly from [74] and [75] and from 
http: / / cfa-www.harvard.edu/ ps/lists/Supernovae.html. 



[72,73]. 
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Figure 5. Limits on the skewness of dark energy arising from the SNIa data when 
flm = 0.3. In the LHS, w = —1 and RHS w = —1.2. The contours correspond to 
68.3, 90, and 95.4 percent confidence Umits. The crosses mark the location of the best-fit 
models (there are two in each figure since, because we marginalize over the direction of 
the coordinates, there is a symmetry of projections about the ^ + 7 = line). 



improved error bars can be used to rule out tliis possibility, and to distinguish whether the 
possible statistical anisotropy was already there at last scattering or whether it is due to 
dark energy. 

5. Inhomogeneous Cosmology 

5.1. Covariant Perturbations 

In the previous sections we have investigated the behaviour of an anisotropic universe 
which is homogeneous. Clearly this cannot be a realistic description of the universe, 
though it suffices as a good approximation when considering the expansion dynamics 
relevant to e.g. the supernovae luminosity distances. When considering e.g. the formation 
of galaxies and other structures in the universe, one has to allow perturbations about 
the homogeneous cosmological solution. Perko et al [76] pioneered the study of galaxy 
formation in anisotropic cosmologies. They discovered that in an anisotropic background, 
the gravitational waves can couple with the density modes. One of their results was that the 
shear can enhance the growth rate of matter perturbations. Later Bianchi I perturbations 
have been also considered within the coordinate approach to cosmological perturbations. 
The perturbation equations have been written quite explicitly in a gauge-ready form 
[77], including then more general sources. However, most studies have considered dust 
cosmologies [78,79]. In view on inflationary models, limiting matter sources to a slowly 
rolling scalar field is of course also of a particular interest. Indeed, very recently there 
have been studies of perturbations in anisotropic inflation which allow an anisotropically 
expanding background [80,81]. Considering late accelerating cosmologies, perturbation 
theory for elastic dark energy has been developed [82]. Yet this has only been done 
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Figure 6. Constraints arising from the SNIa data in the axisymmetric case. Inside the 
darker isosurfaces, the fit is as good as in the ACDM model, < 158. Inside the Hghter 
shaded isosurfaces, one has Ax^ < 8.02. In the left hand side figm-e, 6 — 'y and thus w is 
the equation of state along the symmetry axis. In the right hand side figure, 7 = and 
thus w+S is the equation of state along the symmetry axis. One notes that larger skewness 
6 would typically be compatible with the SNIa data for " supcrnegative" equations of state 
w < —1 and large matter densities. 

in the statistically isotropic framework, though the models would allow an anisotropic 
generalization [24]. In this paper we work within the covariant formalism. Therefore our 
notation follows most closely Dunsby's study of Bianchi I perturbations [83] in the covariant 
approach, which we slightly generalize from the point of view of possible non-isotropizing 
cosmologies, thus taking into account imperfect sources [84]. Maartens and Tsagas have 
considered imperfect sources in the form of magnetic fields in the covariant approach in 
detail [85, 86]. We will omit most of the details of the covariant formalism and refer to, for 
instance the recent review of large-scale structure in relativistic cosmology discussing also 
more general spacetimes than the FLRW model [87]. 

5.1.1. Definitions. In a homogeneous background universe it is then convenient to 
consider variables which are defined as spatial gradients of some physically interesting 
quantities. These should then, by definition, count as perturbations. Since these gradients 
vanish in the background, the variables are by construction gauge-invariant. The basic 
variables we employ are defined as 

Da = -('^V,p, Ya = ^^^VaP, Za = i^^'^VJ, Qa = ^^'^V^^.fe (82) 
P 
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The first two of these are fiuid quantities related to the energy density and pressure. is 
the gradient of the expansion and Vt^ is the rotational perturbation. 

In addition to these, we define some quantities related to the imperfect properties of 
the fluid: 

n„ = h\^\.^,, (83) 

5a = /^^(^c.O^,, (84) 

S„ = £(2) V„a. (85) 

It is also useful to define two scalars, the divergences of the acceleration and of the heat 
transfer, and two vectors, the spatial gradients of these as 

>l = aV, ^„ = (^)V„A g = g%, ga = ^'V„Q. (86) 

Prom these definitions it follows that 

^''V^g, = Q- aaq"", ^''^V\ab = - K^a'^n,,. (87) 

Note that our conventions are slightly different from the Tsagas et al [87], and when 
comparing the results the relations (87) arc useful. 

For a fluid with an equation of state p = p{p, s), where s is the entropy, we decompose 
the spatial gradient of the pressure as 



£Ya = cipDa + pEa, K = - (^^ j V„s, (88) 

where the adiabatic sound speed is, at the zeroth order, related to the equation of state 
w = p/p by 

w = (w-cl){l+w){e + 0. (89) 
We have defined an auxiliary variable ^ to characterize the anisotropy, 

e - (90) 

Using this variable will considerably simplify the following equations. 

Note that the dot corresponds to the derivative along the flow as stated in Eq.(2). It 
corresponds to the derivative with respect to cosmic time t only for scalars. 

5.1.2. Nonlinear Equations. The Raychaudhuri equation is 

9 + \9^ + 2{a^-uj^)-A + l{p + 3p) = 0. (91) 

By projecting the covariant divergence of the stress-energy along and orthogonal to 
m", one obtains, respectively, the two conservation equations 

p+{p + p)e+{l + w)p^ + Q^O, (92) 
(p + p)a„ + y„ + n„ + Q„ = 0, (93) 

where we have defined 

Qa^qa + aJ'qb + ^eqa. (94) 
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The shear evolution equation can be written as 

2 1 

We will be able to drop away the RHS of this equation at the hnearized level. The <> 
means now the traceless symmetric par. Eab is the electric part of the Weyl tensor. This 
far our equations have been exact and thus very general but not very solvable in practice. 



5.1.3. Linearized Evolution of Density. We now linearize the equations, supposing an 
anisotropic but homogeneous background. We assume that the heat flux and vorticity 
vanish at the zero order. All the possible spatial gradients vanish in the background, since 
it is homogeneous. However, we allow anisotropy both in the metric, a ah-, and in the fluid 
TTaft. Note that in this subsection we do not yet specify the type anisotropy, and thus the 
following equations apply in any Bianchi type spacetime. 
We find that the overdensity evolves as 

Da = wdDa - (1 + w)Za - a\D, + (1 + w)^ {Da - lUa) + ^ (^^ " S„) + ^ {OQa - Qa) 

where w = p/ p. The evolution equation for the gradient of the expansion scalar is 

Za = ~pDa - ^eZa - a\Z, - Uo'aa + M„ - 4aE, + + (96) 
From the continuity equation (93) we see that A^, as defined in (86), may be given as 

Aa = -^T^ (('^V^n + (3)V„(^)V^(^) V%,e + (^)V.(3)V^Q.) - 2cle^a 

p(l + ly) ^ 
where one can show that 

(^V.(^)V''Q, = ((^) V,Q)- + 2^(^) V,Q + o'}^^V,Q + 2a^^(^) V,(^V,ge. (97) 

The evolution equation for the density gradients follows then by differentiating the 
first order equation and using all the previous formulas. We find 



Da + 



(- + c^-2«;)e-(l + 2He 



Da + 2(7„^ A 



(98) 



D„ 



+ 



- 2w + cl)9 - {1 + 2w)i 



[l + ci + w)i +[6i + -e{2 + cli-l + 2,cl - 3w) + ^w)]^ 



+ (cf - w)i' 



Da 



E„ 



- w^Ea + aa'Eh) + 2(1 + wyjQa - 4a(l + w)Ea = Si^^ + 5f ). 
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The right hand side of Eq.(98) vanishes if we consider a case that the imperfect 
perturbations can be neglected. The sources are due to the heat flux perturbations as 



(9) 



= - Qa - [(2 + Cl)e + (1 + cDe] Qa - O^Q^ 



+ 



(3w - 2)a2 + p + ^ + (3 + + w)ei + (1 + c2)e 



Qa 



The sources from perturbations of the anisotropic stress are the following 



(n) 



+ 



(3)v„(3)vK3)v%,,-5„ 



+ + (1 + w)^ 



+ 



3(1 - + cf )p + 'ic\a' + ^ + (3 + Ic'M + (1 + <W 



(99) 



(100) 



We cannot derive the general evohition equations for these terms. Their behaviour depends 
completely on the nature of the fluid. In fact this is true also for the other quantities, but 
then we have a simpler and conventional way of parameterizing the unspecified properties 
in terms of c^. 



b.l.Jf.. Linearized Evolution of Shear. In addition to the unspecified matter part, the shear 
(7 and its gradient appear in our evolution equations. So we should find a way to solve 
for them too. Prom the shear evolution equations (95) we easily find that now 

1 

^=a 2 

Prom the definition of E^, and the commutation rules for the spatial and time derivatives, 
one gets that Sa = ^OEa + ^^^^Vafi - eh^^{a.cU'');b + ih'^aih''d^;b);cU''. Inserting the result 
(101) and linearizing then gives 



'''^«f>(3)v7 ^ ^^ab„ f 2/5^2 ^ „ab ( j-. 

-a ^'V(aab)--a aafaf,--Oa --a [Eab--7Tab 



(101) 



(7^ 



^a^V^a^, + lea' + ^a'^ (e,, - ^tt,, 



(102) 



_L "'"(3)V7 -'-^cd(3)r7 „ ^ „cd„ „f o„2 ^ ^cd 



a 



E. 



cd 



f.d 



+ a'Jlc + (.eaaa = 0. 

However, it is also possible to solve the shear away without evolving it in time but relating 
it to the other variables. This is because the time-space components of the Einstein field 
equations yield three additional constraint equations which we have not exploited this far. 
These constraints can be put in the form 



(103) 



As a consistency check, one might take the time derivative of this equation and see whether 
it is compatible with the evolution equation for the shear we just wrote down in (102). It 
turns out that this constraint is indeed preserved in time. 
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We can use this information to eliminate from the matter evolution equation (98). 
In terms of the quantities appearing in that equation, we have 



We note that does of course not vanish even in the general statistically isotropic case. 
5.2. Structure Formation 

In the case of magnetic fields it has been learned that neglecting the background anisotropy 
and resorting to an almost-FLRW treatment can be a good approximation on large scales 
when the magnetic field is weak, but the shear may introduce small-scale effects which 
are accurately captured only by the Bianchi I treatment [86]. As we consider a general 
anisotropic fluid not necessarily reducible to magnetic field plus perfect matter, it is thus 
interesting to study under which conditions we can safely regard the FLRW description 
as accurate, and under which conditions the background anisotropics should be taken into 
account in order not to miss possible new physical effects. 

In the following, we write down the general equations of the last sections in three cases. 
To close the system of equations, one must specify the anisotropic properties of the fluid. 
In subsection (5.2.2) we assume and study the case that the anisotropy is small. Then, it 
turns out, the almost-FLRW approach is valid, with however some direction-dependence of 
the perturbations. In subsections (5.2.3) and (5.2.4) we do not assume that the anisotropics 
are negligibly small at the background order, but make other simple assumptions about 
the properties of perturbations to fully specify the system. In these cases, the shear effects 
could not be described by the almost-FLRW treatment. 

5.2.1. Scalar Equations. Since the matter density can be described as a scalar field, it is 
usual to concentrate on the scalar modes of perturbations when studying linear evolution 
in cosmology. To make contact with these studies, it is useful to extract the scalar part of 
the more general equations in the previous section. A conventional procedure for this is 
based on a local decomposition [88] of the spatial gradients. 



As an example, we may act on the density gradient of matter field, D^, with (^)V", to 
get the spherically symmetric part of the scalar perturbation. Explicitly, we define this 
variable as 




+ - \{9+ ma - <^:^i> +\ma- Qo) - ^.^q. 

p 6 6 




(105) 



^^'^VDa = A. 



(106) 
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In the presence of shear there arise other relevant scalars quantifying the properties of 
matter distribution. We define them as follows 

^a''^^^^VaD^ = A'^, (107) 

iTT^^^^^VaDb = A'^, (108) 

£(7"V^^^V„A = A'^'. (109) 

It then follows that 

f^^V'ba =A + A^, (110) 

e^^^W'^Da = A + 2A^ + + A^' - A'^, (111) 

£(3)vV„^A = A" + + A'"' - A'^. (112) 



In a completely analogous way, one may derive scalar variables from other gradients 
appearing in our considerations. The scalar equations one then derives by acting on the 
general equations with ^(^^V". It is straightforward to pick up the new coefficients of the 
scalar variables by using the definitions (106-109) and the identities (110-112). 

In the most general case, the system is very complicated. In the following, we 
will comment the physical implications under three different approximations and the 
motivations for such approximations: I) The case that the anisotropy is small and can be 
treated as perturbation, II) the case that the from of anisotropy perturbation is determined 
by it's background form III) the case that the the anisotropy is homogeneous also at the 
perturbative level. 



5.2.2. Case I: Small Anisotropy. We make the simpUfying assumption that the anisotropy 
is small. This is justified by the constraints on the late-time anisotropies arising from 
the quadrupole pattern in the CMB and the directional dependence of the luminosity- 
distance. Thus we regard both the metric a and the fluid tt as first order quantities in 
the perturbative equations. Then we are allowed to drop a considerable number of terms 
from the equations, but the largest effects from anisotropy are still taken into account since 
some first order terms remain. Note that would be absent or of a different form in the 
FLRW case. Thus this seems as a very reasonable starting point for the first analysis of 
the structure formation with an imperfect dark energy scenario. 
With this assumption the Eq.(98) reduces to 



;i + 4^-^^^-3c^)p + c^(^V 



Da (113) 



= - w^^^v'^Ea - 2(1 + w)clena + s-f ) 

where the imperfect fiuid terms are given as 

^Cn) ^ i ((3)V„(3) v^(3) V%,c + etla) +3i{l-w + cl) n„. (114) 

Directly from the metric, the anisotropic terms would contribute only quadratically. Prom 
the fluid they contribute also at the flrst order affecting then directly the evolution of the 
density perturbation as the source term (114). 
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By extracting the scalar part of this system as described in subsection (5.2.1), one may 
verify that the result is equivalent to the evolution equation for generalized matter when 
imperfect fluid terms are taken into account [89] . One thus arrives at the conclusion that the 
law governing the evolution of overdensities is exactly the same as in an isotropic case when 
we have neglected the anisotropy terms at the background order. However, one should be 
careful to not deduce that this means that the total evolution is isotropic in such a case. The 
system Eq.(113) presents three equations for the time evolution of three independent spatial 
gradients of the density, and in the general case they could have both different solutions 
and different forms. In particular, this could be a consequence of an intrinsic anisotropic 
stress source, Ila, which could be zero for a given direction a but not for the others. 
Our formula Eq.(113) thus provides a tool to study fully consistently and covariantly the 
fundamentally anisotropic properties of matter inhomogeneities in cosmology, which may 
be hidden in the standard description of the overdensity scalar perturbation. One may 
witness the anisotropy by monitoring the solutions of the individual equations of the system 
(113), while summing their gradients yields the scalar quantity which has a straightforward 
interpretation as the overdensity in the standard FLRW cosmology, which at small scales 
does not depend on the gauge choice. 

This also shows that, even when a Bianchi universe is close to the FLRW model, 
one cannot consider the anisotropics in the perturbed Bianchi universe as a sum of 
FLRW perturbations and the small anisotropics of the exact Bianchi model. When an 
imperfect source is present, its perturbations (114) can break explicitly the statistical 
isotropy already at the linear level. These fluid perturbations of course couple to all other 
perturbations in the metric and in matter. These become then statistically anisotropic as 
well. Therefore the fluctuations of the model evolve differently than in a FLRW model, 
in particular the resulting anisotropy pattern is not a simple superposition of standard 
isotropic perturbations and a Bianchi type perturbation about the FLRW. However, this 
kind of approximation could be reasonably good when only perfect fluid sources are present. 
This is because the shear in the metric typically enters the equations only trough o"^, 
i.e. quadratically. In such a case, one could expect that some effects arc given, to first 
perturbative order, as the sum of ordinary FLRW perturbation and the Bianchi pattern. 

We conclude that when the anisotropy is small, perturbations in each direction behave 
like in different isotropic universes, each obeying their Eq.(113). 

5.2.3. Case II: "Adiabatic" Anisotropy. Clearly, to perform actual calculations one has 
to specify the properties of the fluids under consideration. We would still like to keep the 
discussion general, and therefore we apply a parametrization which could approximate a 
wide range of different models at a certain limit. We then choose the w, 5, 7 parametrization 
as it seems a natural and straightforward generalization of the concept of equation of state 
in the FLRW universe. In the Bianchi I background, one may fully specify the properties 
of a fluid with energy density p by giving the three pressures in different directions, or 
equivalently, the three equations of state for the fluid. In the following, we assume this 
holds to the perturbed order. In principle the inhomogeneities in the stresses would bring 
additional degrees of freedom, but our parametrization assumes they are described by the 
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homogeneous equations of state. This may be interpreted as a generahzed adiabaticity 
condition, since in absence of entropy, the pressure perturbation evolution is determined 
by the homogeneous background evolution. 

Recall the decomposition of the fluid content. In Eq.(15) we described the energy 
momentum tensor for the anisotropic component. It may be written as 

rj^^ = diag (-1, - - 7, + 25 - 7, - + 27) px- 

We are now in a frame having — [1, 0, 0, 0). Therefore, the components in the stress- 
energy tensor (6) could now be identifled as 

P = px, P = wxPx, qa = 0, (115) 
n\ = diag (0, -5 - 7, 25 - 7, -5 + 27) px- (116) 

From the observation that vanishes one could justify to use the equation (98) with the 
heat flux sources 5'(^) turned to zero. If the anisotropy is to be small, the condition is 
clearly that 7 and 6 are small. One flnds that the gradients of the anisotropic stresses arc 
now 

^a = j [0, -(5 + ^)D,, {2S - ^)Dy, {-5 + 2^)D,] , (117) 
and the sum term appearing in the source term Eq.(lOO) is then 

£(3)v„(3)v^(3)v'^7r5e = - (5 + 7)^')VJ3)V-D, (118) 
P 

+ {26 - -ff^Vj^^V^Dy -{6- 2^f^Vj^^V'D,. 

The latter source term thus couples to Dt-tcrms with a ^ b into the evolution equation for 
Da, (98). Since the metric shear = ^ {R^ + — RS) 9^ clearly the small anisotropy 
limit corresponds to -R, S* ^ 1 . 

Let us be more precise about this limit. One observes that if the coefficients 6 and 7 
are supposed to be perturbatively small, all contributions to the perturbations equations 
at the linear order vanish. This is because we consider that the constant coefficient 
of the homogeneous anisotropy gives the relation between the density gradient and the 
(anisotropic) pressure gradient as well as the relation between the background density 
and the homogeneous (anisotropic) pressure. Since the density gradient is already a 
perturbation, the anisotropic pressure gradients (115) should be of the second order in 
perturbations. Similarly with the metric shear: if S and R are regarded as first order, then 
all the terms involving a in the evolution equations arc higher than first order. Therefore 
the anisotropic terms would decouple from all perturbations equations and the effect of 
anisotropy would indeed be negligible. Note that the homogeneous anisotropy and the 
gradient of the anisotropy could be both of the linear order, though it is not the case 
with the present parametrization. Now the simplest nontrivial corrections to the standard 
FLRW treatment appear if one considers the squares of the anisotropy variables as first 
order perturbations. 
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In general, if the anisotropy terms satisfy such conditions, the system of three equations 
for the gradients of the matter aggregations can be written as 

Da+ + - 2«;^ Da + 2(7„^A - Q + 4«; - ^w'' - 3c^^ pD^ (119) 
- Q + 2^ - c^) OaJ-D, - {clDa + + 2(1 + wylOn^ 



= 4(t(1 + w)S„ + - 
P 



(3) v„(3) V^(^) V%,e - 5„ - ( ^ + ^ 1 



+ ^rr„ + 3 (1 - + c^) plla 

Now, in addition to the imperfect source stresses, the couphng of the metric shear will 
generate anisotropies in the perturbations. 

We conclude that there exists a generahzation of the adiabatic condition to Bianchi I 
universes which then uniquely determine the perturbative stresses of the fluid. 



5.2.4.- Case III: Smooth Anisotropy. Yet one more physically motivated approximation 
is to consider the anisotropy field as smooth. A dark energy -like component in usual 
cases is considered to be rather smooth, especially at small scales. In fact no evidence for 
clumping of dark energy has been found. Slowly rolling light fields with minimal couplings 
tend to have negligible perturbations. For the cosmological constant this becomes an exact 
covariant statement: there are no cosmological fluctuations in the usual A-term. When the 
pressure varies in direction and with time, its spatial gradients may be only approximately 
zero. Basically, one then considers the effects of the anisotropic background expansion 
on the matter inhomogeneities. To do this consistently one should couple not only the 
metric a but also the fluid tt anisotropy to the equations. Then we can, however, neglect 
the possible entropy in the perturbations by a direct analogy with the usual A term. One 
notes that the entropy couples also to the large scale perturbations through the anisotropic 
terms, whereas in an FLRW universe one may neg lect E in the hmit (3)v2 0. 

Indeed, for the model of subsection (5.2.2), one could consider a very tidy set of 
perturbation equations where the isotropic perturbations of the A-term identically vanish, 
and the condition ^ = allows to drop a large number of terms. For the variant of the 
model, discussed in subsection 3.2.2, the isotropic pressure of the A-term is not constant 
and therefore cannot be smooth. This introduces entropy between the A-term and matter. 
One would then consider the evolution of the total density gradient Da, as that is simply 
related to the matter density gradient D^'' by Da""^ = (1 — U)Da- The anisotropic stress 
and the metric shear at the background level would then be 

7r"6 = rfm^[0,-l,2,-l]5pA (120) 

a\ = Uiag[Q,l,-2,l]HR, (121) 

giving then cr^ = {HRy/3. Note that tt^ is constant (when S is), and thus its fluctuations 
Ua vanish. However, similar conclusion does not hold for the cr-terms: the metric shear 
has fluctuations and generally is nonzero. Therefore in the model with the constant 
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anisotropic pressure field generates effects on the first order fluctuations. In fact this 
is already reflected in the fact that the total sound speed squared of the system (61) 
can become nonzero when 5 is nonzero. Because there is entropy present, the possible 
negativity of the sound speed does not necessarily lead to instabilities. These effects could 
be used to constrain further the scenario proposed in subsection 3.2.2. 

We conclude that either the isotropic or the anisotropic pressure of an imperfect fluid 
always fluctuates. 

6. Conclusions and Discussion 

We investigated the possibility that the accelerated expansion of our universe, both during 
the early inflation and the present dark energy domination, is driven by a component with 
an anisotropic equation of state. The motivation for this investigation comes from the 
frequent appearance vector fields and possible Lorentz violations in fundamental theories, 
from the need to the test the robustness of the basic assumptions of cosmology, and from 
the hints of statistical anisotropy of our universe that several observations seem to suggest. 
From the examination of the background dynamics we learned that 

• There exist (anisotropic) scaling solutions even without coupling between the 
components. 

• An anisotropic generalization of the A-term has nontrivial dynamics given if its 
isotropic pressure is time- varying. 

• The quadropole constraint on the anisotropy may be avoided, and the present SNIa 
data allows large anisotropy. 

Prom the analysis of system of perturbations the following was found 

• Even in the limit of weak anisotropy, the linearized perturbations will be direction- 
dependent. 

• The linearized anisotropy is the sum of FLRW perturbations and Bianchi I background 
only for perfect-fluid cosmology. 

• Either in the limit of vanishing entropy or in the limit of homogeneous field the almost 
FLRW treatment is not vahd. 

Let us elaborate on these and other results. 

We started by presenting a phenomenological description of a such an anisotropy in 
terms of the skewness parameters 5 and 7, and derived the basic equations describing 
the model as a dynamical system, and obtained the generic asymptotic evolution of the 
universe in some simple cases. Several scaling solutions were found both within the matter 
dominated epoch as well as during an accelerated expansion phase (with the summary 
of fixed points presented in the Table 2). This seems to be interesting, since while in 
the FLRW universe it has been proven difficult to address the coincidence problem by 
finding a model entering from a matter-dominated scaling solution to an accelerating 
scaling solution, allowing for the presence of three expansion rates now opens up the 
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possibility of describing a universe entering from a perfect fluid dominated scaling to 
an anisotropically accelerating scaling era, which might eventually help to understand 
the coincidence problem. In that framework, it is easier to find models with natural 
parameter values avoiding the fine-tuning problems of cosmological constant or scalar field 
quintessence [20, 48, 90]. 

We also contemplated upon the possibility of an anisotropic generalization of the 
cosmological constant constant term. Interestingly, we found that a model with a constant 
density and constant anisotropic pressures (p = po, 8 = 5o, 7 — 7o) is indeed a consistent 
covariant modification of the Einstein's cosmological term which could generate a viable 
expansion history (with small enough 5q and 70) that exhibits nontrivial anisotropic 
properties. We mentioned noncommutative properties of spacetime as a possible origin of 
anisotropics in a cosmological A-term. One could consider different motivations by writing 
the vacuum energy as an integral of the zero-point energies over the k up to some cut-off 
(summed over every possible field) . If there were any isotropy-breaking physical principle 
which allowed only wavevectors with specific directions, the vacuum energy would acquire 
nonzero 5 and 7. More specifically, if the cut-off depended not only on the magnitude k of 
the wavevector k, but also on its direction (i.e. some directions would be excluded from 
the integral or contribute with less weight to the vacuum energy) this could explain both 
the smallncss of the observed cosmological constant and the apparent anisotropic features 
in the cosmological data. 

The amount of possible anisotropy, as described by the skewness parameters, can 
be constrained by the CMB and SNIa observations. We found that the present CMB 
data constraints can be much tighter than the constraints from the anisotropics in the 
luminosity-distance -redshift relationship of SNIa. However, there are classes of models 
where the CMB bounds from the background quadrupole anisotropy can be avoided 
altogether while there is significant anisotropy at late times that could be observed in 
the late SNIa distributions and which could leave statistically anisotropic imprints at the 
large-angle CMB. Wc point out, however, that in more general models, e.g. with time- 
varying 5 and 7, one could allow more anisotropy. It is in principle possible for an arbitrarily 
anisotropic expansion to escape detection from CMB (considering only the effects from the 
background), as long as the expansion rates evolve in such a way that — ey — In 
other words, the (background) quadrupole vanishes, if each scale factor has expanded - no 
matter how anisotropically - the same amount since the last scattering. 

We also considered the inhomogeneous cosmology of these models. The covariant 
formalism was employed to derive the general evolution equations for the perturbations in 
a universe with background anisotropics, i.e. imperfect fiuid terms and shear of the metric. 
Our formalism and equations thus provide a tool to study fully consistently and covariantly 
the fundamentally anisotropic properties of matter inhomogeneities in cosmology, which 
may be hidden in the standard description of the overdensity scalar perturbation. One 
may witness the anisotropy by monitoring the solutions of the individual equations of 
the system (113), while summing their gradients yields the scalar quantity which has a 
straightforward interpretation as the overdensity in the standard FLRW cosmology, which 
at small scales does not depend on the gauge choice. 
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The implications of anisotropic sources for the formation of structures were discussed 
in several limiting cases. In particular, in the cases examined, the deviations from isotropy 
are small enough that one may treat all the imperfect terms as first order perturbations. We 
noticed that this gives a result which is not obtainable by a straightforward generalization of 
the usual FLRW equations. If one considers the evolution equation of matter perturbations 
(in its closed form (113)), one notices that while the metric shear would enter only at the 
second order, the imperfect fluid terms are present at the first order level. Thus, even 
when the Bianchi universe would be close to the FLRW model, one cannot consider the 
anisotropics in the perturbed Bianchi universe as a sum of FLRW perturbations and the 
small anisotropics of the exact Bianchi model. This is seen from the fact that when an 
imperfect source is present, its perturbations (114) can break explicitly the statistical 
isotropy already at the linear level. In general, the perturbation structure of these models 
allows a rich array of possibilities. The anisotropic properties seem to usually extend to 
the smaller scales also, but this feature is model-dependent and the constraints it implies 
must be studied on a case-by-case basis. 

The present SNIa data allows anisotropic acceleration, but future SNAP data could 
set things straight about the skewness of dark energy and so of its nature. Such possibility 
would open a completely new window not only on the nature of the CMB anomalies but 
also into high energy physics models beyond the usual isotropic candidates of inflation or 
dark energy models such as the inflaton, the quintessence scalar fields or even an isotropic 
cosmological constant. Hence, an anisotropic equation of state may not only successfully 
unify the early inflation with the present days acceleration, but might also be the culprit 
for both the apparent problem of cosmic coincidence as well as the large-angle anomalies 
in the CMB. 

The detailed signatures of the anisotropics, either originating from the early or the late 
accelerating epoch of the universe, to the full sky maps of the CMB should be determined: 
could these signatures rule some or all of the models out, or could they in fact be identified 
with the observed anomalies? Investigations along these lines are ongoing. 

7. Appendix: Pedestrian notation for Bianchi I 

In this Appendix we review our system within the coordinate formalism. Our main 
developments in the bulk of the paper follow the covariant approach for the sake of 
homogeneity. However, for some considerations the full covariant machinery is unnecessary, 
and analysis is simpler and perhaps more transparent by specifying the coordinate system 
and writing the equations in terms of the metric components. 

The Bianchi I metric as generalization of the flat FLRW metric is 

ds^ = -dt^ + a^{t)dx^ + b^{t)dy^ + c^{t)dz^. (122) 

Due to the presence of three scale factors a, h and c, one has to introduce also three 
expansion rates. In principle all these could be different, and in the limiting case that all 
of them are equal, one of course recovers the FLRW case. Note that we use the roman 
letters to label the covariant indices, like in the expression for the energy momentum tensor 
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(6), whereas we use the Greek letters to label the indices of tensors when we refer to the 
components in the specific choice of basis (122). The perfect fluid energy momentum tensor 
is then simply 

'^(m)v = diag{-l, Wm, Wm, Wm)Pm, (123) 

In the coordinate system we have chosen, the energy-momentum tensor must be diagonal, 
since the nondiagonal Einstein tensor vanishes identically. However, the pressure in each 
spatial direction could be different. We then parametrize the anisotropy of the pressure by 
the two skewness parameters 

35 = {py - Px)/p, 37 = (p^ - p^)/p. (124) 

The energy-momentum tensor of this fluid is then 

T"^^ diag{-l,w,w + SS,w + 3-f)p. (125) 

The Einstein equations would not allow non-diagonal components for the stresses with the 
metric (122). 

With the matter content (123) and (125), the three nontrivial spatial components of 
the field equations 



are given as 



b c be 

- + - + - = -8nG [wmpm + wp] , (126) 
b c be 

(i c cic 

- + - + - = -SttG \WraPm + (^^ + 35) p\ , (127) 

a c ac 

Qj h Qjb 

- + - + — = -SttG [Wmpm + (W + 37) p] , (128) 

a b ab 
and the time component is 

ab be cd ^ ^ , , /-,.^^\ 

- + - + — ^8nG[pm + p]. 129 
ab be ea 

Overdots denote derivatives with respect to time t. 

The continuity equations are then given by the divergence of the energy-momentum 

tensors. We let the two components also interact, and thus allow nonzero divergence for 

the individual components. 

Pm + 3H{1+ Wm) Pm = QHp, (130) 

and 

p + H[3{l + w) + 6{3-2R + S)+-f{3 + R- 2S)] p = -QHp, (131) 

where Q describes the coupling. If there are no interactions, one finds that the densities 
scale as 

Pm ~ (a6c)-^-"'-, p ~ (a6c)-^-'"fc=^*c^^. 
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The time variable we will use is an average e-folding time defined as follows: 

x=-\og{ahc). (132) 
3 

One notes that H — x,\n analogy with the FLRW case, where the e-folding time reduces 
to X — logo. Derivative of a function / with respect to x will be denoted by prime, unless 
/ is explicitly specified as a function f{y) of some other variable y. 

In the following we will rewrite this system in a more convenient form for analytical 
and numerical study. It is then useful to follow the notation introduced in Ref. [49] by 
expressing the mean expansion rate as an average Hubble rate H 

1 I d b c\ , ^ 

and the differences of the expansion rates as the Hubble-normalized shear R and S 

Assuming that all the expansion rates are positive, one has always — 3 < i?, 5' < 3. This is 
true also if either i? = or 5" = 0. We find then a generahzed Friedmann equation of the 
form 

^ - 3 i-i(R2 + s^-Rsy 

Our dynamical variables are then the density fraction U, and the two shear anisotropics 
R and S. The evolution equations for these are then 

U' ^U(U -l)[^(3 + R-2S) + 5(3-2R + S) + 3(w- w^)] - UQ, (136) 
S' ^ ^{9 - R^ + RS - S^) {S [U {S + J + w - Wm) +Wm-1]- 67C/} ,(137) 

R' ^^(Q-R^ + RS-S^){R[U{6 + 'y + w- w^) + w^-l]- 66U} .(138) 
Acknowledgment s 

DFM acknowledges support from the A. Humboldt Foundation. We thank the anonymous 
referee for critical reading of the manuscript and several suggestions which have helped us 
to improve the paper. 



References 

[1] H. K. Eriksen, F. K. Hansen, A. J. Banday, K. M. Gorski, and P. B. Lilje. Asymmetries in the cmb 

anisotropy field. Astrophys. J., 605:14-20, 2004. 
[2] G. Hinshaw et al. Three-year wilkinson microwave anisotropy probe (wmap) observations: 

Temperature analysis, astro-ph/0603451, 2006. 
[3] Kate Land and Joao Magueijo. The axis of evil. Phys. Rev. Lett., 95:071301, 2005. 
[4] Alcksandar Rakic and Dominik J. Schwarz. Correlating anomalies of the microwave sky: The good, 

the evil and the axis, astro-ph/0703266, 2007. 



CONTENTS 



40 



[5] Joao Magueijo and Rafael D. Sorkin. Occam's razor meets wmap. Mon. Not. Roy. Astron. Soc. Lett., 
377:L39-L43, 2007. 

[6] Kate Land and Joao Magueijo. The axis of evil revisited, astro-ph/0611518, 2006. 

[7] Craig J. Copi, D. Huterer, D. J. Schwarz, and G. D. Starkman. On the large-angle anomalies of the 

microwave sky. Mon. Not. Roy. Astron. Soc, 367:79-102, 2006. 
[8] Amir Hajian and Tarun Souradeep. Testing global isotropy of three-year wilkinson microwave 

anisotropy probe (wmap) data: Temperature analysis. Phys. Rev., D74:123521, 2006. 
[9] Christopher Gordon, Wayne Hu, Dragan Huterer, and Tom Crawford. Spontaneous isotropy breaking: 

A mechanism for cmb multipole alignments. Phys. Rev., D72:103002, 2005. 
[10] Cristian Armendariz-Picon. Footprints of statistical anisotropics. JCAP, 0603:002, 2006. 
[11] A. E. Gumrukcuoglu, Carlo R. Contaldi, and Marco Peloso. Cmb anomalies from relic anisotropy. 

astro-ph/0608405, 2006. 

[12] Lotty Ackerman, Sean M. Carroll, and Mark B. Wise. Imprints of a primordial preferred direction 

on the microwave background. Phys. Rev., D75:083502, 2007. 
[13] Stephen H. S. Alexander. Is cosmic parity violation responsible for the anomalies in the wmap data? 

hep-th/0601034, 2006. 

[14] C. Armendariz-Picon. Creating statistically anisotropic and inhomogeneous perturbations. 

arXiv:0705.1167 [astro-ph], 2007. 
[15] Christian G. Bohmer and David Fonseca Mota. Cmb anisotropics and inflation from non-standard 

spinors. arXw:0710.2003 [astro-ph], 2007. 
[16] Angelica de Oliveira-Costa and Max Tcgmark. Cmb multipole measurements in the presence of 

foregrounds. Phys. Rev., D74:023005, 2006. 
[17] Kaiki Taro Inoue and Joseph Silk. Local voids as the origin of large- angle cosmic microwave 

background anomalies. Astrophys. J., 648:23-30, 2006. 
[18] Kari Enqvist and Tcppo Mattsson. The effect of inhomogeneous expansion on the supernova 

observations. JCAP, 0702:019, 2007. 
[19] Christian Armendariz-Picon. Could dark energy be vector-like? JCAP, 0407:007, 2004. 
[20] C. G. Boehmer and T. Harko. Dark energy as a massive vector field. Eur. Phys. J., C50:423-429, 

2007. 

[21] Maxim Libanov, Valery Rubakov, Eleftherios Papantonopoulos, M. Sami, and Shinji Tsujikawa. Uv 
stable, lorentz-violating dark energy with transient phantom era. JCAP, 0708:010, 2007. 

[22] Tomi S. Koivisto and David F. Mota. Vector field models of inflation and dark energy. 0805.4229, 
2008. 

[23] Tomi Koivisto and David F. Mota. Dark energy anisotropic stress and large scale structure formation. 

Phys. Rev., D73:083502, 2006. 
[24] Richard A. Battye and Adam Moss. Anisotropic perturbations due to dark energy. Phys. Rev., 

D74:041301, 2006. 

[25] P. R. Pereira, M. F. A. Da Silva, and R. Chan. Anisotropic self-similar cosmological model with dark 

energy. Int. J. Mod. Phys., D15:991-999, 2006. 
[26] Kishore N. Ananda and Marco Bruni. Cosmo-dynamics and dark energy with a quadratic 

eos: Anisotropic models, large-scale perturbations and cosmological singularities. Phys. Rev., 

D74:023524, 2006. 

[27] Bijan Saha. Anisotropic cosmological models with perfect fluid and dark energy revisited. Int. J. 

Theor. Phys., 45:952-964, 2006. 
[28] Bijan Saha. Anisotropic cosmological models with perfect fluid and dark energy. Chin. J. Phys., 

43:1035-1043, 2005. 

[29] Valeria Pettorino, C. Baccigalupi, and F. Perrotta. Scaling solutions in scalar-tensor cosmologies. 
JCAP, 0512:003, 2005. 

[30] Valeria Pettorino, Carlo Baccigalupi, and Gianpiero Mangano. Extended quintessence with an 

exponential coupling. JCAP, 0501:014, 2005. 
[31] Tomi Koivisto. The matter power spectrum in f{R) gravity. Phys. Rev., D73:083517, 2006. 
[32] Tomi Koivisto and Hannu Kurki-Suonio. Cosmological perturbations in the palatini formulation of 



CONTENTS 



41 



modified gravity. Class. Quant. Grav., 23:2355-2369, 2006. 
[33] Timothy Clifton, David F. Mota, and John D. Barrow. Inhomogeneous gravity. Mon. Not. Roy. 

Astron. Soc, 358:601, 2005. 
[34] Constantinos Skordis, D. F. Mota, P. G. Ferreira, and C. Bochm. Large scale structure in bekenstein's 

theory of relativistic mond. Phys. Rev. Lett., 96:011301, 2006. 
[35] Morad Amarzguioui, O. Elgaroy, D. F. Mota, and T. Multamaki. Cosmological constraints on f{R) 

gravity theories within the Palatini approach. Astron. Astrophys., 454:707-714, 2006. 
[36] David F. Mota and Douglas J. Shaw. Strongly coupled chameleon fields: New horizons in scalar field 

theory. Phys. Rev. Lett., 97:151102, 2006. 
[37] David F. Mota and Douglas J. Shaw. Evading equivalence principle violations, astrophysical and 

cosmological constraints in scalar field theories with a strong coupling to matter. Phys. Rev., 

D75:063501, 2007. 

[38] Tomi Koivisto. Growth of perturbations in dark matter coupled with quintessence. Phys. Rev., 
D72:043516, 2005. 

[39] Luca Amendola. Coupled quintessence. Phys. Rev., D62:043511, 2000. 

[40] Robert Caldwell, Asantha Cooray, and Alessandro Melchiorri. Constraints on a new post-general 

relativity cosmological parameter, astro-ph/0703375, 2007. 
[41] Kiyotomo Ichiki and Tomo Takahashi. Constraints on generalized dark energy from recent 

observations, astro-ph/0703549, 2007. 
[42] D. F. Mota, J. R. Kristiansen, T. Koivisto, and N. E. Groeneboom. Constraining dark energy 

anisotropic stress. arXiv:0708.0830 [astro-ph], 2007. 
[43] Anthony W. Brookficld, C. van dc Bruck, D. F. Mota, and D. Tocchini-Valcntini. Cosmology of 

mass-varying neutrinos driven by quintessence: Theory and observations. Phys. Rev., D73:083515, 

2006. 

[44] Luca Amendola, Martin Kunz, and Domenico Sapone. Measuring the dark side (with weak lensing). 

arXiv:0704.2421 [astro-ph], 2007. 
[45] Carlo Schimd, Jean-Philippe Uzan, and Alain Riazuelo. Weak lensing in scalar-tensor theories of 

gravity. Phys. Rev., D71:083512, 2005. 
[46] D. F. Mota and C. van dc Bruck. On the spherical collapse model in dark energy cosmologies. Astron. 

Astrophys., 421:71-81, 2004. 
[47] Nelson J. Nunes and D. F. Mota. Structure Formation in Inhomogeneous Dark Energy Models. Mon. 

Not. Roy. Astron. Soc, 368:751-758, 2006. 
[48] Tomi Koivisto and David F. Mota. Accelerating cosmologies with an anisotropic equation of state. 

arXiv:0707.0279 [astro-ph], 2007. 
[49] John D. Barrow. Cosmological limits on slightly skew stresses. Phys. Rev., D55:7451-7460, 1997. 
[50] John D. Barrow, Yoshida Jin, and Kei-ichi Maeda. Cosmological co-evolution of yang-mills fields and 

perfect fiuids. Phys. Rev., D72:103512, 2005. 
[51] Luigi Campanelli, P. Cea, and L. Tedesco. Ellipsoidal universe can solve the cmb quadrupole problem. 

Phys. Rev. Lett, 97:131302, 2006. 
[52] Jose Beltran Jimenez and Antonio L. Maroto. Cosmology with moving dark energy and the cmb 

quadrupole. astro-ph/0703483, 2007. 
[53] Davi C. Rodrigues. Anisotropic cosmological constant and the cmb quadrupole anomaly. 

arXiv:0708.1168 [astro-ph], 2007. 
[54] Michael J. Longo. Does the universe have a handedness, astro-ph/0703325, 2007. 
[55] Michael J. Longo. Is the cosmic axis of evil due to a large-scale magnetic field, astro-ph/0703694, 

2007. 

[56] John D. Barrow and Sigbjorn Hervik. Anisotropically inflating universes. Phys. Rev., D73:023007, 
2006. 

[57] John D. Barrow and Sigbjorn Hervik. On the evolution of universes in quadratic theories of gravity. 

Phys. Rev., D74:124017, 2006. 
[58] C. P. Burgess, Richard Easther, Anupam Mazumdar, David F. Mota, and Tuomas Multamaki. 
Multiple inflation, cosmic string networks and the string landscape. JHEP, 05:067, 2005. 



CONTENTS 



42 



[59] Damien A. Easson, Ruth Gregory, David F. Mota, Gianmassimo Tasinato, and I. Zavala. Spinflation. 

JCAP, 0802:010, 2008. 
[60] George F. R. Ellis and Henk van Elst. Cosmological models, gr-qc/9812046, 1998. 
[61] Roman V. Buniy. Arjun Berera, and Thomas W. Kephart. Asymmetric inflation: Exact solutions. 

Phys. Rev., D73:063529, 2006. 
[62] Bijan Saha. Interacting spinor and scalar fields in bianchi type-i universe filled with viscous fluid: 

Exact and numerical solutions, gr-qc/0703124, 2007. 
[63] Bijan Saha. Bianchi type i universe with viscous fluid. Mod. Phys. Lett, A20:2127-2144, 2005. 
[64] I. Brevik and S. D. Odintsov. On the cardy-verlinde entropy formula in viscous cosmology. Phys. 

Rev., D65:067302, 2002. 

[65] I. Brevik and O. Gorbunova. Dark energy and viscous cosmology. Gen. Rel. Grav., 37:2039-2045, 
2005. 

[66] Tver Brevik. Crossing of the w = -1 barrier in two-fluid viscous modified gravity. Gen. Rel. Grav., 

38:1317-1328, 2006. 

[67] Emma ,1. King and Peter Coles. Dynamics of a magnetised bianchi i universe with vacuum energy. 

Class. Quant. Grav., 24:2061-2072, 2007. 
[68] Bali Raj, Pareek Umesh Kumar, and Pradhan Anirudh. Bianchi type i massive string magnetized 

barotropic perfect fluid cosmological model in general relativity. arXiv:0704.0753 [gr-qc], 2007. 
[69] Andrew Pontzen and Anthony Challinor. Bianchi model cmb polarization and its implications for 

cmb anomalies. arXiv:0706.2075 [astro-ph], 2007. 
[70] Shin'ichiro Ando and Marc Kamionkowski. Nonlinear evolution of anisotropic cosmological power. 

arXiv:0711.U779 [astro-ph], 2007. 
[71] Cyril Pitrou, Thiago S. Pereira, and Jean-Philippe Uzan. Predictions from an anisotropic inflationary 

era. arXw:0801.3596 [astro-ph], 2008. 
[72] Camillc Bonvin, Ruth Durrer, and M. Alice Gasparini. Fluctuations of the luminosity distance. Phys. 

Rev., D73:023523, 2006. 

[73] Dominik J. Schwarz and Bastian Weinhorst. (an)isotropy of the bubble diagram: comparing 

hemispheres. arXiv:0706.0165 [astro-ph], 2007. 
[74] Adam G. Riess et al. New bubble space telescope discoveries of type ia supernovae at 2: > 1: Narrowing 

constraints on the early behavior of dark energy, astro-ph/0611572, 2006. 
[75] Pierre Astier et al. The supernova legacy survey: Measurement of omegam, omegaiambda and w 

from the first year data set. Astron. Astrophys., 447:31-48, 2006. 
[76] T. E. Perko, R. A. Matzner, and L. C. Shepley. Galaxy formation in anisotropic cosmologies. Phys. 

Rev., D6:969-983, 1972. 

[77] H. Noh and J. C. Hwang. Perturbations of an anisotropic space-time: Formulation. Phys. Rev., 

052:1970-1987, 1995. 

[78] P. G. Miedema and W. A. van Leeuwen. Perturbations in bianchi i universes. Class. Quant. Grav., 
9:S183-S185, 1992. 

[79] E. Dimastrogiovanni, W. Fischler, and S. Paban. Perturbation Growth in Anisotropic Cosmologies. 

arXiv:0803.2490 [hep-th], 2008. 
[80] A. E. Gumrukcuoglu, Carlo R. Contaldi, and Marco Peloso. Inflationary perturbations in anisotropic 

backgrounds and their imprint on the cmb. arXiv:0707.4179 [astro-ph], 2007. 
[81] Thiago S. Pereira, Cyril Pitrou, and Jean-Philippe Uzan. Theory of cosmological perturbations in an 

anisotropic universe. arXiv:0707.0736 [astro-ph], 2007. 
[82] Richard A. Battye and Adam Moss. Cosmological perturbations in elastic dark energy models. Phys. 

Rev., D76:023005, 2007. 

[83] Peter K. S. Dunsby. Gauge invariant perturbations of anisotropic cosmological models. Phys. Rev., 
D48:3562-3576, 1993. 

[84] Peter K. S. Dunsby. Gauge invariant perturbations in multicomponent fluid cosmologies. Glass. 

Quant. Grav., 8:1785 1806, 1991. 
[85] Christos Tsagas and Roy Maartens. Magnetized cosmological perturbations. P/jys. i?et;., D61:083519, 
2000. 



CONTENTS 



43 



[86] Christos G. Tsagas and Roy Maartens. Cosmological perturbations on a magnetised Bianchi I 

background. Class. Quant. Grav., 17:2215-2242, 2000. 
[87] Christos G. Tsagas, Anthony ChaUinor, and Roy Maartens. Relativistic cosmology and large-scale 

structure. arXiv:0705.4397 [astro-ph], 2007. 
[88] G. F. R. Ellis, M. Bruni, and J. Hwang. Density gradient - vorticity relation in perfect fluid robertson- 

walker perturbations. Phys. Rev., D42:1035-1046, 1990. 
[89] Tomi Koivisto. Viable palatini-f(r) cosmologies with generalized dark matter. Phys. iiew., D76:043527, 

2007. 

[90] Jose Beltran Jimenez and Antonio L. Maroto. A cosmic vector for dark energy. arXiv:0801.1486 
[astro-ph], 2008. 



